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Abstract 

We consider (4, l)-dimensional branes constructed with two scalar fields <f> and x coupled to 
a Dirac spinor field by means of a general Yukawa coupling. The equation of motion for the 
coefficients of the chiral decomposition of the spinor in curved spacetime leads to a Schrodinger-like 
equation whose solutions allow to obtain the masses of the fermionic modes. The simplest Yukawa 
coupling fy(j)x^ is considered for the Bloch brane model and fermion localization is studied. We 
found resonances for both chiralities and related their appearance to branes with internal structure. 
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I. INTRODUCTION 



Braneworld scenarios have received great attention 
dressing important problems such as gauge hierarchy 



iom the physical community after ad- 
l|, [2I, [3] and the cosmo logical constant 



% |5] . The construction of domain walls embedded in a higher dimensional bulk, however, 



has a more ancient history in the literature 



As is well-known, (4, 1)-D branes can 



be classified as thin and thick ones. Thin branes are constructed after introducing a tension 
te rm in the action, locked by a Dirac de.ta function (see Q and fences therein). The 
mathematical formalism that describes the metric around the brane allows to obtain several 
interesting results in the domain of high energy physics and cosmology. Usually, the issue of 
field localization in such branes is addressed with the help of Dirac delta functions, without 
any clear subjacent dynamics. On the other hand, thick branes are constructed in a dy- 



namical way after introducing one or more scalar fields coupled with gravity ll|, |l2j, |l3j, [14]. 



Thick branes are more natural in the sense that field and gravity localization can be studied 
with the introduction of smooth functions (instead of Dirac ones). Moreover, the thin brane 
solutions are recovered in certain limits Q]- One important issue in (4, 1)-D braneworld 
models is to consider spacetimes with a bulk metric that extends the Randall-Sundrum re- 
sult [l7] in a way that gravity can be localized in the (3, 1)-D slice [lO, 3, Jil 20, 21]. The 
spacetime around the Randall-Sundrum brane is anti-de Sitter (AdS^,), and several thick 
brane models (that present this characteristic asymptotically) have been constructed. 



Braneworld models with one or more scalar fields were constructed in order to re 
gravity localization with brane thickness 22j, |23j. In particular, the Bloch brane model 



26|, 



ate 



27] 



is composed of two scalar fields and is the extension for the known Bloch walls 
in the context of braneworlds. Bloch branes are stable structures that can be characterized 
as branes with interfaces located at the maximum of energy density. The two interfaces 
signals the presence of an internal structure in the brane. The presence of internal structure 
in topological defects was firstly considered in Witten's superconducting cosmic strings 24]. 
There the internal structure is provided by the condensation of a (charged) field over the 
defect. In the Bloch brane case, the appearance of internal structure is concerned with a low 
value of the coupling constant that guides the way the two fields interact with each other. 
For a stronger interaction, a single peak for the energy density characterizes a brane with 
simpler structure. It was also shown 23] that in the presence of gravity the interfaces are 



2 



located more closely in the extra dimension. As such branes are able to localize gravity, one 
interesting question is investigating how fermionic fields may be localized. In this way, there 
arises the expectation that a brane with richer internal structure could provide new results 
when compared to a simple kinky brane. 

lie study of fermion localization on branes 
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36|, 133, |38J, 



is rich and interesting. In order to localize fermions in branes one needs a coupling 
between the spinors and the scalar fields that form the brane. This is a condition also 
present in Jackiw and Rebbi treatment of fermion localization on solitons in flat space 41]. 
In branes, the procedure consists in separating from the full spinor one scalar coefficient 
with dependence only on the extra dimension, leading to a Schrodinger-like equation and 
a probabilistic interpretation. Depending on the model, one can obtain resonant massive 
states. 

In particular, the authors of j^] analyzed the issue of fermion localization on a brane 
constructed with the sine-Gordon potential. That work considered general fermionic Yukawa 
couplings between one scalar field and spinorial fields. It was found that the simplest Yukawa 
coupling ^0^/ allowed left-handed fermions to possess a zero-mode that localizes on the 
brane. The right-handed fermions, on the other hand, present no zero-mode. The large 
massive fermionic modes (for both left and right chirality) are plane waves and thereby are 
not localized on the brane. 

Inspired on the results obtained in 37j], we investigate the issue of fermion localization 
of fermions on a brane constructed from two scalar fields coupled with gravity. In Sec. 
2, we present the known first-order equations for a brane model coupled with scalar fields 
Sec. 3 deals with the fermionic sector of the model, where we extend the results 



1_W( 



from [37|] for two scalar fields <fi and \ with general coupling. Afterwards, we choose the 
simple coupling ^ipx^ ■ In Sec. 4, we consider the Bloch brane model and our main results 
are presented. The Schrodinger-like equation for the model is obtained and numerically 
solved. The solutions correspond to massive fermionic modes. We have found resonances on 
the brane and have related their appearance with the simultaneous increasing of the brane 
internal structure. Our conclusions are presented in Sec. 5. 
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II. SCALAR FIELD AND METRIC EQUATIONS 



Our system is described by the action 



S = I d^xdyy/W\[~R+^d a ^ + ^d aX ^x-V(M], 



(1) 



and the metric 



ds 2 = g a bdx a dx b , 

ds 2 = e 2A r) IJ , 1/ dx* 1 dx u — dy 2 , 



(2) 



where g = det (g a b)- Here a,b — 0,1,2,3,4, and e is the warp factor. We suppose that 
A = A(y),<ft = <f)(y) and \ = x{v)- The extra dimension y is infinite, and we have the 
continuity of the first derivatives of (f)(y), xiv) an d A(y) with respect to y as boundary 
conditions for the thick brane. 

The action given by Eq.(0Q) leads to the following coupled differential equations for the 
scalar fields <p(y), x(y) an d the function A(y) from the warp factor: 



0" + 4A'<p' 



dV^,x) 

dV(<P,x) 
dx 



(3) 



(4) 




(5) 



(6) 



where prime stands for derivative with respect to y. 



With the potential [ill ] . 





the first-order differential equations which also solve the equations of motion are 




(8) 



1 dW 



(9) 



2 d x ' 

A' = --W. 



(10) 
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III. FERMIONIC SECTOR 



Now, we consider a Dirac spinor field coupled with the scalar fields by a general Yukawa 
coupling. The action for this sector is 

S = J d 5 xy/=g{^r M D M y -v^F((f),x)% (11) 

Here we consider the fields <j) and % as background fields given as solutions of the Eqs. (J3j)- (EI) - 
In other words, we neglect the backreaction from the spinor field on the brane solutions. 
We change variable from y to z with 

z = f V A «)^, (12) 
Jo 

in order to rewrite Eq. ([2]) and get a conformally flat metric 

ds 2 = e^^^dx^dx" -dz 2 ). (13) 

This change of variable is usual in problems of gravity localization. Indeed, for a large class 
of models, one can achieve a Schrodinger-like form for the equations for metric fluctuations, 
when decoupled from the scalar fluctuations in a specific gauge. In the treatment of fermion 
localization the same change of variable is used. The equation of motion corresponding to 
the action given by Eq. ( fTTl) is 

[y^ + + 2d z A) - r]e A F(<P, = 0. (14) 

We apply the well known procedure of a general chiral decomposition for the spinor 

(x)a Ln (z) + ^2ip Rn (x)a Rn (z), (15) 

n n 

where the sum can be over discrete bounded modes or over a continuum of modes. In 
this decomposition ipLn(x) and ipR n (x) are, respectively, the left-handed and right-handed 
components of the 4-dimensional spinor field with mass m n . Also, there are two scalars ai n 
and «fl n that depend only on the extra dimension z. In order to turn the notation easier, 
we dropped the index n in which follows. After applying Eq. ( fi~5i) in Eq. ( fi~4"l) . we obtain 
two equations for the scalars and o>r. 

[d z + 2d z A + V e A F{<P, X )}a L {z) = ma R {z), (16) 
[d z + 2d z A - r/e A F(0, X )Wr(z) = -ma L (z). (17) 
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Implementing the change of variables ai = e 2A L and a R = e 2A R 37|, we get Schrodinger- 
like equations for the wavefunctions R m (z) and L m (z) given by 

H L L m (z) = m 2 L m (z), (18) 
H R R m (z) = m 2 R m (z), (19) 

with corresponding Hamiltonians 

"* = (-s +e V)(s + «v)- < 20 > 

The structure of the Hamiltonians guarantees that m is real, the absence of tachyonic modes 
and the possibility of interpreting \R m (z)\ 2 and \L m (z)\ 2 as the probability for finding the 
right and left massive modes in a given coordinate z. The Schrodinger equations can be also 
written explicitly as 

- L" m {z) + V L (z)L m (z) = m 2 L m} (22) 
-Rl{z) + V R {z)R m (z) = m 2 R m , (23) 

with 

V L (z) = e 2A V 2 F(<f>, x? ~ e A V d z F(4>, X ) - (d z A)e A V F(<f>, X ), (24) 
V R (z) = e 2 V^(0, X? + e A V d z F(<f>, X ) + (d z A)e A V F^, X ), (25) 

named as the Schrodinger-like potentials for the fermionic fields. The Hamiltonians from 
Eqs. (1201) and (12 ip can be written as Hi = A^A and H R = AA^ . This shows that Hi and 
H R are conjugated Hamiltonian of supersymmetric quantum mechanics and Vi(z) and V R (z) 
are superpartner potentials. In this way there is a zero mode for Vi(z) and the spectra of 
Vl(z) and V R (z) are interrelated (42, 43]. 

From now on, we consider the simplest Yukawa coupling F(<p, X ) = 4>X an d study the 
occurrence of massive modes for the known Bloch brane model with two scalar fields. 



IV. THE BLOCH BRANE MODEL 



A rich model that leads to branes with internal structure is 

2 



25, 



26 



27] 



- 2a<p X z 



(26) 



y 



FIG. 1: (a) Warp factor e 2A ^ (left) and (b) field (f>(y) (right) for a = 0.05 (thinner trace, larger 
curves), a = 0.10 and a = 0.20 (thicker trace, narrower curves). 



with a being a real parameter (0 < a < 0.5). In this case, 
can be solved analytically, leading to the following results 

(j){y) = tanh(2ay), 

x(y) = 



he first-order equations (I8~f)-(1T0~ 
23J: 



and 



2 ) sech(2aw/), 



A(y)=— (1 - 3a) tanh 2 (2ay) - 21ncosh(2ay) 
9a 



(27) 
(28) 

(29) 



Fig. [T^i shows how the warp factor tends to zero going to the AdS^, limit far from the 
brane. The brane thickness decreases with a, and this can be seen by the plots of the field 
solutions in Figs. [Qd and [2^. 

Now, we turn to the z- variable according to Eq. f[T^]) . For general values of a, the 
expression for A(y) from Eq. f)29p is not suitable to be integrated in a known explicit form 
and numerical methods are necessary. However, for the particular case a = 1/3 we have 



A(y)=-- In cosh(2ay), 



and 



cosh( -y 



(30) 



(31) 



After simple integration, we get an expression for z(y) in an explicit form in terms of an 
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FIG. 2: (a) Field x{u) (left) and (b) function z(y) (right) for a = 0.05 (thinner trace, larger curves), 
a = 0.10 and a = 0.20 (thicker trace, narrower curves). 



hypergeometric function: 



y_ 

\y\ 



9 / , / 2 \ \ 2/3 / 112 / ,/2\\'\ 9 



TlcoshU,)] 2Fl I cosh -_r - T - 



(32) 

Unfortunately, even in this simpler case, we can not obtain the inverse y(z) in an explicit 
form. However, this explicit form for z(y) is already useful to produce sequences of pairs 
of points (z, y) with a constant step in z. Indeed, given a value z the corresponding y 
can be found after numerically solving the implicit equation (Eq. (1321) ). This leads to 
the determination of the functions <p(z) , x{z) , A(z) with good precision. The corresponding 
derivative with respect to z, A' (z), can be found as 

dA = dAdy = dA / 2 \ ~ 2/3 
dz dy dz dy \ 3 / 

with dA/dy stemming from the explicit expression for A(y) (Eq. ( 1301) ). Similarly, we can 
construct the functions (f)'(z), x'( z )- All the former calculations can be used to construct 
the Schrodinger potentials Vl and Vr displayed in Fig. [3] with good precision; this was used 
with the purpose of validating the numerical method used for general values of a. . 

For some values of a, we integrated numerically Eq. f|T2|) obtaining z(y), whose graphic 
is depicted in Fig. [2Jd. With this, we can numerically determine A(z), <p(z), x( z )- m 
order to obtain the derivatives A'(z), (j)'(z), x'{ z )i a numerical procedure was constructed for 
generating sequences with constant step in z. In this way the procedure for general values 
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FIG. 3: Comparing procedures: Schrodinger-like potential (b) Vl(z) (left) and (c) Vr(z) (right) 
for a = 1/3, obtained with explicit function z(y) (points) and with the full numerical procedure 
(line). Note that both procedures give almost indistinguishable results. 
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FIG. 4: Schrodinger-like potentials (a) Vl{z) (left) and (b) Vr(z) (right) for a = 0.05 (thinner 
trace), a = 0.10 and a = 0.20 (thicker trace). 

of a is more involved once we do not have an explicit form for z(y) to guide us as in the 
a = 1/3 previous case. Eqs. ( 124|) and ( 1251) are then used to provide a graphical form for 
the potentials Vl and Vr. We applied this procedure for the a = 1/3 case, whose results 
are depicted in Figs. Ek-b. Such graphics reveal that both procedures described here lead 
to almost the same results. 

Now, we can use the general numerical procedure to investigate the potential Vl for other 
values of a. Some results are exhibited in Fig. H^i. We note from this figure that Vl has 



a volcano-like shape, that goes asymptotically to zero far from the brane. This form of 
potential is usually found in problems involving gravity localization. Here, we can search 
for the existence of a zero mode followed by a continuum of massive modes by means of an 
analogy. Indeed, we note that the potential is negative at the brane location; this guarantees 
the existence of a normalized zero mode after integrating Eq. (lisp for m = 0: 

L ocexp(-7// dz'e A ^()){z')x{z')). (34) 
Jo 

Fig. Hb shows that the potential Vr is always positive at the brane position, which is not 
compatible with the existence of bound fermions with right chirality. However, we note the 
appearance of a hole in the potential that grows for lower values of a. This fact is new and 
our hypothesis is that it can be related, in some way, to the internal structure 22| of the 



brane constructed with two fields. It is remarkable that such behavior of the Vr potential 
was not observed in a previous treatment for one-field models. The appearance of this hole 
in the potential can be responsible for resonances or at least for a light increase in the decay 
rate of massive fermions on the brane. Note also that the hole around z = for Vr is absent 
for larger values of a, in principle prohibiting the existence of resonances for this range of 
parameter. In order to look for resonant effects, we must consider now the massive modes, 
solving numerically the Schrodinger-like equation with purely numeric potentials. 



A. Right-handed fermions 

Firstly we will consider the case of right-handed fermions with the small parameter a = 
0.05. We can vary m? to get the wavefunctions R m (z). Eq. (1231) shows that we need two 
initial conditions. As we are dealing with thick brane with no tension at z = 0, the Israel 
junction conditions applied to the brane give R(z) and R'(z) continuous at z — 0. Note 
that the absence of negative values for the potential Vr prohibits the existence of bound 
states. However, the tunneling process across the barrier (localized at the region between 
the maximum of the potential and the origin) can result in different rates for the leaking 
process depending on the mass of the modes. Here the region of interest corresponds to 
masses satisfying 0.42 < m 2 < 2.44, where the tunneling process is effective (see Fig. Hb - 
thinner trace). 

Figs. [5^-d shows typical even massive fermionic modes before normalization for the region 
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FIG. 5: Non-normalized even massive right-handed fermionic modes for a = 0.05 with (a) m 2 = 1 
(upper left), (b) m? = 2 (upper right), (c) m 2 = 2.1 (lower left) and (d) m 2 = 4 (lower right). 

z > 0. Figs. [5^,-c belongs to the more interesting region 0.42 < m 2 < 2.44 whereas Fig. 
[5jd is for a higher massive mode. The lower massive modes show that there is a transient 
behavior followed by typical plane wave oscillations, characteristic of a free mode. This 
shows that these modes represent massive fermions that certainly will be leaked from the 
brane. Comparing Fig. [5^ with Figs. [5)d and[5fc, we note that the transition region roughly 
corresponds to —20 < z < 20. Note also that for lower values of m 2 the values of the 
non- normalized wavefunctions R(z) depart largely from the initial value 1. More important 
is to estimate how larger is the solution R(z) in the transition region —20 < z < 20 in 
comparison to the amplitude of the plane wave oscillations. Indeed, one would expect that 
broader differences would correspond to greater lifetimes for the massive fermionic states 
near to the brane. This difference assumes a huge amount for m 2 = 2.1, as we can see from 
Fig. Efc. 
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FIG. 6: (a) Normalized squared wavefunction for right-handed fermion on the brane, |i?(0)| 2 , as 
a function of m 2 . Note the high peak characteristic of resonance. For the present case we have 
a = 0.05. (b) Normalized squared wavefunction for right-handed fermion on the brane, |iZ(0)| 2 , as 
a function of m, identifying better the region defining the peak of resonance. 

This lead us to the point of normalization of the massive modes and to the better estimate 
of R m (0). Their importance is that we can rewrite Eq. {23} as ] R R R{z) = m 2 R(z), which 
ensures that we can interpret |i? m (0)| 2 as the probability for finding the massive modes on 
the brane. In particular, large peaks in the distribution of -R(O) as a function of m would 
reveal the existence of resonant states closely related to the existence of a brane with internal 
structure. 

We can consider a normalization procedure for the wavefunctions R m {z) in a box with 
borders located far from the turning point, where the solutions have characteristics of plane 
waves toward the brane. In this way, we considered a box with —100 < z < 100. An 
important information can be obtained from the normalized |i? m (0)| 2 , since from it we can 
compare the relative probability for finding a massive mode on the brane. The result is 
depicted in Fig. EK We note from the figure a huge peak around m 2 = 2.1, characterizing 
the occurrence of a resonance and a long-lived massive fermionic mode on the brane. We 
note from the figure that low massive modes have negligible associated probabilities, whereas 
high massive modes are characterized by a plateau in the probability distribution. This is 
related to the finite size of the box used for normalization. Indeed, for m 2 » Vr we 
can approximate the Schrodinger equation as —R"{z) ~ m 2 R(z), with normalized solution 
R{z) ~ 1 / ^/z m ax cos(mz) , where 2z max is the size of the box used for normalization. For 
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FIG. 7: (a) Normalized squared wavefunction for right-handed fermion on the brane, |i?(0)| 2 , as a 
function of m 2 , for a = 0.10 and —100 < z < 100. (b) Normalized squared wavefunction for right- 
handed fermion on the brane, |i?(0)| 2 , as a function of m, identifying better the region defining the 
peak of resonance. 

our choice z max = 100, we have |i?(0)| 2 = l/(z max ) = 0.01, corresponding to the plateau 
observed in Fig. Ek 

Since in our model the extra dimension is infinite, one may wonder about the influence 
of z max on the normalization of the spectral density. Indeed, another choice of z max would 
lead to another plateau in the probability distribution. The physical information, however, 
is contained in the value of the resonance peak, which does not depends on z max since it is 
chosen sufficiently large. 

Fig. [6] shows the peak of resonance with the horizontal scale now taken as m. We can 
estimate the life-time r for the observed resonance from the width at half maximum of the 
peak appearing in the figure. Indeed, the width in mass V = Am of the peak is related to the 
life-time. In this way, in its own reference frame, the fermion disappears toward the extra 



dimension with time scale r ~ T" 1 [44]. The peak was located at m — 1.451(1) with the 
maximum probability around 0.109. The width at half maximum is Am = 0.0083, resulting 
in a lifetime r a = ro.05 = 120. These results are in arbitrary units, and their significance will 
be better understood after comparison with the results for other values of the parameter a. 

Now, we repeat our analysis for other values of a to investigate the effects of the coupling 
parameter a between the fields <fi and \ for the existence of resonances and their lifetimes. 

For a = 0.10 the results for |i?(0)| 2 , as a function of m 2 , are depicted in Fig. [7k Note 
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a=0.10 



a=0.05 



a=0.20 




FIG. 8: Normalized squared wavefunction for right-handed fermion on the brane, |-R(0)| 2 , as a 
function of m 2 . We used —100 < z < 100, and the resonant peaks correspond to a = 0.05 (third 
narrow and smaller peak), a = 0.10 (higher peak). Also it is showed the broader peak for a = 0.20. 

from the figure that we also have one resonant peak as before. Now, however, the larger 
value of a results in a smaller value of the mass for the resonant state. Also, the resonant 
peak is more pronounced. Fig. [7)d shows that the peak was located at m = 0.968 with the 
maximum probability around 0.202. The width at half maximum is Am = 0.02, resulting 
in a lifetime r a = r 01 = 50. Comparing this result with r a = r 05 = 120, we conclude that 
smaller values of a result in resonances with larger lifetimes. 

We compare both cases in Fig. [HI where we display the resonances for a = 0.05 (second 
smaller peak) and a = 0.10 (first larger peak). We know that smaller values of a correspond 



to branes with richer internal structure 



23]. Then the observed thinner peaks for smaller 



parameters a means resonances more pronounced and a trapping mechanism more effective. 
This is connected with the larger amount of matter now forming the brane due to the 
coupling of the scalar field <fi to the other field \- 

The reduction of the lifetime for resonances with larger values of a poses the question if 
there are no resonances at all on the brane for values of a above a certain threshold a*. We 
noted the decreasing of the associated mass for the resonances with larger values of a. Also, 
it was noted a qualitative change on the potential Vr above a certain value a*. For example, 
for a = 0.20, Fig. [9^ shows that now the potential assumes the form of a volcano. One 
characteristic of this type of potential is the occurrence of a zero- mode when the central part, 
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FIG. 9: (a) Schrodinger-like potential Vr(z) for a = 0.20. Note the appearance of a region with 
negative potential, absent for the lower values of a considered previously, (b) Normalized squared 
wavefunction for right-handed fermion on the brane, |i?(0)| 2 , as a function of m. Note now the 
presence of a broader peak with Am > m. 



around the brane, assumes negative values. This is the potential observed in Fig. for 
left-fermions. Here, however, the central part of the potential is positive and the zero-mode 
is absent. For a = 0.20 there were observed no resonances on the brane. In fact, Fig. [9b 
depicts the probability |i?(0)| 2 as a function of the mass of excitations. Here we can see that 
the width at half maximum Am is larger that the position of the peak m = 0.95(1). This 
characterizes a state with a lifetime r = l/Am too small to result in any physical effect, 
characterizing the absence of a resonant state. 

We can relate these results with the behavior of the energy density of the brane without 



fermions as studied in Ref. 



231 ] . The energy density T^(y) as a function of the coordinate 



y and for a fixed parameter a is given by 



T a (y)=e 



2A{r) 



fd<p(y)V (dx(y)X 
[2[\ dy ) V dy ) 



v(<f>(y),x(y)) 



(35) 



Fig. [10] shows that the structure of the brane changes when the parameter is a = 0.17. 
Indeed, we note that for 0.17 < a < 0.5 the energy density characterizes a defect with the 
energy density centered around a central peak, whereas for a < 0.17 the appearance of two 
symmetric peaks signals the occurrence of a brane with internal structure (more details in 



Ref. 



231]). Now, we connect this known result with our findings of resonance peaks for 



small values of a and with the absence of such resonances for larger values of a. There 
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FIG. 10: (a) Plots of the matter energy density T ^ 05 (y) (thinner trace), T ft 10 (y),T ft 17 (y) and 
T 00 2 °(y) (thicker trace). 

appears to be a connection between the occurrence of metastable states on the brane and 
the occurrence of two peaks in the energy density characterizing these branes. From this we 
conclude that branes with internal structure favor the appearance of resonance states for 
right-handed fermions. Also, for branes without internal structure, we were not able to find 
resonance states. 

A further point to remark is that we have considered analytical fittings for the Schrddin- 
ger potentials Vr(z) in terms of finite series with oscillating terms. For fittings that agreed 
visually with our numerical points, the Runge-Kutta method was applied and the normalized 
massive modes were constructed. The results for the resonances were roughly the same, 
corroborating our conclusions. 

B. Left-handed fermions 

For left-handed fermions, as we already noted, there is a normalized zero mode. We 
repeated the analysis from the previous section and also found resonances. However, for 
left-handed fermions there are differences with respect to the position of the resonance peak, 
in comparison to which was obtained for right-handed ones. For instance, for a = 0.05, we 
found a resonance peak around m = 1.59, with Am = 0.07, corresponding to a life-time 
t = 14. When comparing with corresponding results for right-handed fermions (m = 1.45 
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and t = 120), we see that the left-handed resonance is the more massive and the less 
stable. Analysis for the resonance peak for a = 0.10 showed that, similarly to the observed 
for right-handed fermions, an increasing a also thickens the resonance peak for left-handed 
fermions. 



C. Correspondence between the spectra and realization of Dirac fermions 

The lack of correspondence between the spectra of left- and right-handed fermions appears 
to be in contradiction to the fact that Vl and Vr are superpartner potentials. The reason 
for this is that the parity-odd modes do not couple to the four dimensional braneworld at 
z = 0, because for such modes we have -R(O) = and L(0) = 0. This was also noted in the 



case of gravity localization in another model 45j . After we submitted this work, we knowed 
about Ref. [46[, where a study of odd parity modes is done. According to ref. [46|, we set 
L(0) = and L'(0) = 5 for parity odd wavefunctions as our starting point for the application 
of the numerical Numerov method. We use the scheme with step h in z corresponding to 
the points z$ described by 

L(0) = 0, (36) 

L(h) = hL'(0), (37) 

T( , 2[1 + 5Fte-i]Lte-i) - [1 - F(z^ 2 )]L( Zi - 2 ) 

L(z>) - f — ^ , (38) 



where 



F{z) = ~h*<y L {z)-m*). (39) 



The condition 1/(0) = 5 is arbitrary and will be fixed after the normalization process. 

Thick branes are extended objects along the extra dimension that allow us to interpret, 
after normalizing L(z), the probability for finding the massive modes on the brane (not 
necessarily at z = 0) as f_ b z dz\L m (z)\ 2 . The parameter Zb is chosen in order to allow the 
influence of the odd modes in small regions around z — 0. We can consider a normalization 
procedure for the wavefunctions L m (z) in a box with borders located at z = ±z max , far 
from the turning point, where the solutions have characteristics of plane waves far from the 
brane. After normalization of L(z), the probability for finding the massive modes on the 
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FIG. 11: Plots of P brane for z max = 100, z b = 10, a = 0.05. Curves for a) V R (left) and b) V L 
(right) constructed with even (thicker trace) and odd (thinner gray trace) wavefunctions. 



brane is given by 46j 



fl b Zb dz\L m (z) 



Pbrane ~ dz\L m (z)\*' (40) 

J —Zrnax 1 v ' 1 

We followed the same procedure for right-handed fermions and depicted in Fig. [11] the 
results for Pbrane in the particular case of a = 0.05. The values of m for the resonance 
peaks in the figure corresponding to the parity even wavefuntions agree with the obtained 
previously in the plots of the normalized |i?(0)| 2 (Fig. Eb) and |L(0)| 2 . The figure also shows 
that with this procedure the resonance peaks for left- and right-quirality agree with good 
precision and the formation of Dirac fermions is realized. Larger values of the parameter 
Zb do not change the results for the position of the resonance peaks, but small values of Zb 
are more efficient to identify the peaks. One example of such procedure in a model with 
three thin branes was done in Ref. (30], were one of the wavefunctions (right or left-handed) 
is symmetric and the other corresponding to the same mass is antisymmetric with respect 
to the center of the orbifold. This is a consequence of demanding the action to be even 
under the Z2 orbifold symmetry 47J. In this way it is shown the existence of metastable 
massive Dirac modes trapped on the brane. In our case of one thick brane immersed in a 
bulk with infinite extra dimension, if we allow measurements only at z = 0, the evidence of 
Dirac fermions is lost. In the more realistic scenario, one must consider the finite size of the 
brane and a nonzero overlap between the parity odd wavefunctions with the SM fields , 
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which would be sensitive to to the extra dimension in the vicinity of the brane, justifying 
the procedure followed in this section. 



D. Dimensional analysis and phenomenology 



After restoring the dimensional parameters in the action ([TJ), it reads as 



S = I dx 4 dy\/\g 



' -R+lg<*d a <j>d b <f>+lg ab d a xd b x-V(<f>,x) 



(41) 



In a five-dimensional space-time, the mass dimension for the fields and coupling constants 
in the action are: [0] = [x] = M 3 / 2 , [V ((j),x)} = M 5 , [dab] = 1, [R] = M 2 and the five- 
dimensional gravitational constant has dimension \G^ 5 '\ = M~ 3 . 
The potential of the scalar fields (J7J) is written as 

where W is the superpotential with mass dimension 4. The coupling in the last term 
guarantees the possibility of the model to have some compatibility with the expected LHC 
phenomenology. 

In that way, the Bloch brane is defined by the following superpotential 

2/x,o2ci/ t c,o 

W = 2kZ<P - — 3 - — X 2 (43) 
with [Z] = M 3//2 , [k] = M and [a] — 1. It yields the following solutions for the scalar fields 

n \ 1/2 

= Z tanh (2any) , X = Z I 2 j cosh -1 (2any) , (44) 

where Z, the kink amplitude, is given in terms of the five-dimensional gravitational constant 
as Z~ 2 = G^ 5 ' and, the parameter k~ 1 is related with the brane width. Thus, we have only 
two free parameters in the model k and G^ . 

From the explicit form of the potential in term of the scalars fields 

V (0, X) = ^ (Z 2 - 2 - a X 2 ) 2 + - ^|^0 2 [Z 2 - ^ 2 - a X 2 ) , (45) 

we get the masses for the fields and x as being 

[U r- 

rricj, = \ —k , m x = v2an, (46) 
V <j 
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therefore, the parameter k is explicitly related with the field masses. It offers an opportunity 
for k to be the energy scale of the model. 

The action for massless Dirac fermion is given by 

S f = J dx A dy^\ [ $>T a DJ> - 1V VF (0, X ) *] , (47) 

where 7 is a dimensionless parameter and [^f] = M 2 . As the function F(<ft,x) = 4>X has 
dimension 3, we have [rj] = M~ 2 . 

In order to obtain an expression for the effective coupling constant 77 in terms of our two 
free parameters, we check the linear equations of motion (fl6|) and (ITT)) after the change of 
variables ai = e~ 2A L and or = e~ 2A R: 

[0, + 7»7e A F(0,x)]£ = mR } (48) 

[d z - -fr]e A F (0, x)] R = -mL. (49) 
By performing the following rescaling: 

z = H , d z = X-%, (50) 
with [A] = M" 1 , the Eqs. (ji8"|) - (pjl are rewritten as 

[d^ + 7 e A F^,x)}L = mR (51) 

[d^- ie A F{^x)]R = -mL, (52) 

where 4> = <P/Z, % — x/Z are dimensionless fields. We have set 

A" 1 = r]Z 2 , (53) 

and rh = m\ is a dimensionless constant, this way, the parameter A -1 gives our mass or 
energy scale. 

The choosing of the coupling constant rj is a fine-tuning procedure so that the model 
has the chance of providing a phenomenology in agreement with LHC energy scale. Thus, 
setting 77 = kG^\ the mass scale results 

A" 1 = k, (54) 
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therefore the mass scale is related to the brane width and it is independent of the relation 
between and the four-dimensional G constant. 

From the kink equations (j44j) . we define the effective brane width as 

AL = — !— . (55) 

In order for resonances to appear, our model imposes the condition a < 0.17, this gives 
AL > 1.47/c -1 . By considering k~ 1 TeV, we obtain that the brane width is > 10~ 18 cm. 
It is easy to note that if k > 1 TeV we get better lower limit for the brane width. 

V. CONCLUSIONS 

In this work we considered branes constructed with two scalar fields. We considered a 
simple Yukawa coupling between the two scalars and the spinor field. Depending on the 
amount of the coupling parameter a between the scalars, the energy momentum density of 



the branes can be characterized by a two-peak or by a single peak. After |22J we consider 
this two-peak distribution as characterizing a brane with internal structure. For these more 
complex branes, the coupling parameter a is small. We investigated the occurrence of mas- 
sive modes with both chiralities solving numerically the Schrodinger equation and looking 
for zero-modes and possible resonances. For branes with internal structure, we found left- 
and right-handed resonances together with a zero-mode left-handed solutions. The absence 
of zero-mode solutions for right-handed fermions together with their findings for left-handed 
ones agrees with the well-known fact that massless fermions must be single-handed in a 
brane model 



321 ] . With respect to resonances for left-handed fermions, we can resume our 
conclusions with the following points: i) Larger values of a correspond to larger peaks of res- 
onance, with corresponding smaller lifetimes, ii) Above a threshold value of the parameter a, 
the resonances become too unstable and cease to appear. This coincides with the qualitative 
change of the energy momentum density of the brane. iii) Branes with internal structure are 
more effective in trapping left-fermions. iv) The lower is the parameter a, the higher is the 
mass of the the resonant mode. This means that branes with internal structure tend to trap 
matter with larger mass more efficiently in comparison to branes without internal structure. 
A further investigation on this subject could consider other classes of Yukawa interaction in 
order to look for possible effects on the massive fermionic resonances. 
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